Quantum resonances and decay of a chaotic fractal repeller observed using microwaves 
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The quantum resonances of classically chaotic n-disk geometries were studied experimentally utiliz- 
ing thin 2-D microwave geometries. The experiments yield the frequencies and widths of low-lying 
resonances, which are compared with semiclassical calculations. The longtime or small energy be- 
havior of the wave-vector auto-correlation gives information about the quantum decay rate, which 
is in good agreement with that obtained from classical scattering theory. The intermediate energy 
behavior shows non-universal oscillations determined by periodic orbits. 



The n-disk scattering problem is one whose quantum- 
classical correspondence has received extensive theoreti- 
cal attention [|]J^, because it is a paradigm of an open 
quantum chaos system, much as the pendulum is for inte- 
grable systems. Furthermore it is relevant to physical sit- 
uations in diverse fields, such as the crossroads geometry 
for electron devices Q, unimolecular chemical reactions 
and electromagnetic and acoustic scattering Q. 

The classical scattering function of n-disks on a plane 
is non-differentiable and forms a Cantor set. A cen- 
tral property is the exponential decay of an initial dis- 
tribution of classical particles, hence the name (fractal) 
repeller. For closed quantum chaotic billiards, experi- 
mental and theoretical work is available on eigenvalues 
(which are purely real) and eigenfunctions, and the many 
features of universality and non-universal behavior are 
known B. In contrast, for open quantum systems, the 
eigenvalues are intrinsically complex and their universal 
behavior is a question of great interest [||. Despite ex- 
tensive theoretical treatment there have been almost no 
real experiments on the n-disk geometry which exemplify 
this unique problem in quantum chaos. 

In this paper we present a microwave realization of the 
chaotic n-disk problem. Experiments were carried out 
for n = 1, 2, 3, 4, 6, as well as for large n = 20, the latter 
corresponding to the random Lorentz scatterer. In this 
paper we focus on the case n = 4. The experiments yield 
the frequencies and the widths of the low lying quantum 
resonances of the four-disk repeller. We have also carried 
out semiclassical calculations of the resonances, which are 
shown to reproduce the resonances reasonably well. Our 
experiments enable us to explore the role of symmetry in 
a unique way by studying different irreducible represen- 
tations. The experimental data are used to display the 
signatures of the classical chaos in the transmission spec- 
tra, through measures such as the spectral (wave-vector 
k) auto-correlation function. The small k (long time) be- 
havior of this quantity provides a measure of the quantum 
escape rate, and is shown to be in good agreement with 
the corresponding classical escape rate. For large k (short 
times), the contribution of periodic orbits is observed as 



non-universal oscillations of the auto-correlation. 

The experiments are carried out in thin microwave 
structures consisting of two highly conducting Cu plates 
spaced rf ~ 6 mm and about 55 x 55 cm, in area. Discs 
and bars also made of Cu and of thickness d are placed 
between the plates and in contact with them. In order 
to simulate an infinite system microwave absorber ma- 
terial was sandwiched between the plates at the edges. 
Microwaves were coupled in and out using loops termi- 
nating coaxial lines which were inserted in the vicinity of 
the scatterers. All measurements were carried out using 
an IIP8510B vector network analyzer which measured 
the complex transmission (S'21) and reflection (^n) S- 
parameters of the coax -I- scatterer system. It is crucial 
to ensure that there is no spurious background scatter- 
ing due to the finite size of the system. This was verified 
carefully as well as that the effects of the coupling probes 
were minimal and did not affect the results. 

In this essentially 2-D geometry. Maxwell's equa- 
tion for the experimental system is identical with the 
Schrodinger time-independent wave equation (V^ -I- 
k^)"^ = with ^ = Ez the z-component of the mi- 
crowave electric field. This correspondence is exact for all 
frequencies fc < c/2d — 25GiJz.(Note that k = iirf/c, 
where c is the speed of light). It is this mapping which 
enables us to study the quantum properties of the n-disk 
system. For all metallic objects in the 2-D space between 
the plates, Dirichlet boundary conditions apply inside the 
metal. Thus for the 4-disk geometry, ^ = inside the 4 
disks. 

The transmission function S2i{f) which we measure 
is the response of the system to a delta-function exci- 
tation at point fi probed at a different point r2, and 
is determined by the wavefunction \1/ at the probe loca- 
tions riand ^2- In our experiments the coax lines act 
as tunneling point contacts, and hence it can be shown 
that S'2i(/) = A{f)G{ri,f2,f) is just the two-point 
Green's function G(ri,r2,/), scaled by a slowly vary- 
ing function A{f) of frequency / which represents the 
impedance characteristics of the coax lines and probes. 
Because we ensure that the coupling to the leads is very 



weak, any shifts due to the leads are neghgible (< 10~^ 
of the resonance frequencies and widths) 1^ . 

Note that the measurements of \S2i{f)\'^ are equiv- 
alent to a two-probe measurement of the conductance 
Q, since the conductance g cc J2 l^nmP, the sum of the 
transmission probabilities tnm for all input/output chan- 
nels {n,m). In the present experiment, the input/output 
leads are ^-function point contacts, and hence in — 
5'2iand g oc |5'2i(/)p. Thus our experiments enable us 
to make comparison with theories P,0-ll| originally de- 
veloped for electronic micro-structures [l2ln3]. 

The transmission function |<S'2ip for a square 4-disk 
system with disk radius a = 2cm and separation R = 
8 cm is shown in Figjl|. The data shown is the sum of 
traces taken for 3 locations of the probes, in order to 
avoid accidental zeros of the wavefunction at the probe 
locations. Several resonances are clearly seen. The res- 
onance frequencies / and widths A/ directly yield the 
complex wave- vector quantum eigenvalues of the 4-disk 
repeller. In our previous experiments on closed cavi- 
ties such as Sinai billiard geometries jlj], the widths of 
the eigen-resonances were due to dissipation in the metal 
walls. Here however the resonance widths are due to de- 
cay of the wavefunction into infinite space. These exper- 
iments thus enable us to explore quantum decay in open 
systems. It must be emphasized that the dissipation in 
the walls is entirely negligible in the present experiments. 




FIG. 1. Transmission function 1521 P for a 4-disk system 
with R = 8cm and a = 2cm. The dashed line is the 
semi-classical calculation. (Inset) the experimental configu- 
ration and the various reduced representations of the 4-disk 
geometries studied in this work. The separation distance R 
is shown as a dashed line. The x's mark typical locations of 
the coupling probes. 

The experimental approach was validated by measure- 
ments (not shown) on an (integrable) 2-disk system, for 
which we have measured both in the full space with two 
actual disks {Ai and A2 representations) and in the half- 
space using a reflecting mirror {A2 representation). Here 
the spectrum consists of a 1-parameter family of reso- 
nances due to a single periodic orbit. The experimental 
data are in very good agreement with semiclassical cal- 



culations (for preliminary measurements, see |15[|). 

For the chaotic 4-disk case, we have also carried out 
semiclassical calculations using methods in the literature 
[0,p|jl9 pfl] . In the semiclassical theory, the resonances of 
the system are given by the poles of the Ruelle ( function 
with j running from to cx3 p7| 



C(i/2)+,Hfc)=n[i-(-i)^^^'"''/A(^ 



l/2)+J^ 



(1) 



where k is the wave vector, Ip is the length of the periodic 
orbit p , Lp the number of collisions of the periodic orbit 
with the disks, and Ap is the eigenvalue of the stability 
matrix. 

For the 4-disk system with C^y symmetry, the semi- 
classical calculations using the cycle expansions |^ were 
carried out up to period 3 including 14 periodic orbits. 
Because of the C41, symmetry, there are five representa- 
tions of the Ruelle C function |2^, Ai, A2, Bi, B2, E, 
where the last one is two-dimensional. The poles of the 
first Ruelle ( function with j = are calculated, since 
they contribute to the resonances with the longest life- 
times. The energy of the particle is £' = k'^{here we can 
take m = 1/2, h = 1), and is also complex. 

In Fig.n^, we show the results of the semi-classical cal- 
culations, as a superposition of Lorentzians S\'^^i{k) = 
J2i Ci7i/((fc — Si)^ + 7?)i where Si and 7^ are the calcu- 
lated real and imaginary parts of resonances, respectively. 
The coupling constants Ci were chosen to fit the data - as 
explained above they are specific to the probe locations. 
Good agreement is found for the resonance frequencies of 
the sharp resonances. The sharp resonances with small 
imaginary parts and hence high Q are easily recognized 
while the resonances with large imaginary parts are not 
easy to distinguish, although all resonances contribute to 
the transmission function in Fig.l. 

The approximate semiclassical theory used provides a 
fair prediction of the resonance frequencies even for the 
low-lying resonances. The accuracy is within 5%. But 
the agreement for the widths is not as good as that for 
the resonance frequencies. The lower the frequency, the 
greater the discrepancy between the calculated and ex- 
perimental widths. The discrepancy may improve if more 
orbits of higher period are included. Another source of 
the discrepancy is intrinsic in the semi-classical theory 
because of the large correction of the stationary phase 
approximation pl[ |. Comparison with exact quantum 
calculations |[l7|,p| would be desirable. 

One of the noteworthy features of the experiments is 
the ability to vary geometry. We exploited this by ex- 
ploring the role of the symmetry in the 4-disk geometry. 
Wc performed the experiment in four different setups (see 
inset to Fig.ll]). These correspond to four different ways 
of probing the phase space: the full space in which all 
five representations are included, half {A2,B2,E), one 
fourth {Bi,B2), and one eighth of the space (i?2)- Ap- 



proximately 17 configurations were studied and analyzed 
in details. The systematic trends are consistent with ex- 
pectations. Here we discuss some of the principal features 
of the results, further extensive details will be published 
in a longer publication 0. 

We now turn to another analysis of the data in terms of 
the spectral auto-correlation function, which was calcu- 
lated as C{k) =< \S2i{k - {K/2))\^\S2iik + {k/2))\^ >k. 
The average is carried out over a band of wave vector 
centered at certain value fco and of width Afc. The win- 
dow function [^ chosen was f{x) = (1 — |a;|/\/6)/\/6 for 
|a;| < a/6 and /(x) = Ofor |a;| > ^6 with a; = (fc-fco)/Afc. 
The trace 5*21 used is the average of several traces col- 
lected at different probe locations with the same geo- 
metrical configuration of the disks to avoid missing reso- 
nances due to the accidental coincidence of either probe 
with a node of the wavefunction. The autocorrelation is 
the average of that for several fco . A plot of the autocorre- 
lation for the l/8th configuration of the 4-disk geometry 
is shown in Fig.0. 
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FIG. 2. Wave- vector autocorrelation C{k) of the 4-disk sys- 
tem with R — 20cm and a — 5cm. Data are shown for l/8th 
configuration of the 4-disk geometry corresponding to the B2 
representation. The correlation is calculated with interval 
Ak = 2cm~^ . The different sets represent different values of 
the central wave- vector fco. The bold line is a Lorentzian with 
7q„ = 0.075cm"\ 

Since |52i(fc)P=EzC.7^/((fc-5.)2+7'), we have H] 
C{k) = t:J2,^^ c,Cjij, + -/j)/iiK- {si- Sj)y + {-/i + -ijf). 
In the case that there are no overlapping resonances, 
|si — Sjl >> (7i + 7j), the small k behavior of the auto- 
correlation is C(k) ~ TrJ2i 2cf7i/(K^ -I- 47^). 

According to semiclassical theory, the above sum can 
be replaced by a single Lorentzian ||2^,|3|, p5| , p2[ 



C(/t) = C(0) 



1 



1 + («/7)2 



(2) 



In semi-classical theory, 7 = 70, the classical escape rate 
with the velocity scaled to 1. Thus one can interpret 
the width of the autocorrelation as an average lifetime of 
resonances pq,p^. The above equation was used to fit 
the spectral autocorrelation for small n and thus obtain 



the value of the experimental escape rate 7gm, as shown 
in Fig.| 

In the classical scattering theory, the classical escape 
rate 70 is the simple pole of the classical Ruelle C, function 
p6[ . We use the results of ref. ||I^ in Fig. 3 for compar- 
ison with the experimental data. Asymptotically, when 
R is large, 70 « \n{2^/2R/a)/{\/2R). Another relevant 
quantity is the abscissa of absolute convergence Sc for 
Eq. (1) which can also be estimated from the Ruelle ( 
function with the classical cycle weights tp replaced by 
the corresponding semiclassical ones. Sc serves as a lower 
bound of the escape rate |2^, and is also shown in Fig. 3. 
(Note that this is negative for R/a < 4.5). 
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FIG. 3. Experimental escape rate 7,™ scaled to radius 
a= \ versus Rja. Data are shown for different reduced config- 
urations of the 4-disk geometry : l/8th space(open squares), 
1/2 space (open circles), 1/4 space (filled squares), full space 
(triangles). The classical escape rate (solid line) is calculated 
from the first three periodic orbits in the fundamental do- 
main. The abcissa of convergence Sc of Eq.(l) is shown as a 
dot-dashed line, and represents a lower bound on the quantum 
escape rate. 

Good agreement of the escape rate is obtained between 
that of the classical theory 70 and that from the experi- 
mental data 7qm- This is shown in Fig. 3, where we com- 
pare the experimental escape rates 7^™ with the classical 
escape rate 70 for several values of Rja. Note that in 
Fig. 3, data are included for 17 configurations of the dif- 
ferent reduced (1/8, 1/4, 1/2 and full space) representa- 
tions of the 4-disk geometry shown in Fig.l. The radius 
of the disks used was a = hem for the 1/8 space, and 
a ~ 2cm for the others. The data for 7g„j are scaled 
to radius a = \. In Fig. 3 we have shown experimentally 
that the small k, behavior of the spectral auto-correlation 
has a universal behavior in that it is independent of the 
details of the geometry. The good agreement in Fig. 3 
between the measured escape rate 7g,„ and the classical 
escape rate 70 shows that some quantum properties are 
well predicted by semiclassical theory. 

For intermediate k, the semiclassical theory Eq.(2) fails 
because of the presence of the periodic orbits, which lead 
to non-universal behavior. In the case of just one periodic 



orbit, C(k) cx X^^o^t/K^ ~ nAs)"^ + 47^]. For exam- 
ple for the 2-disk problem if just the A2 representation 
is present ||, As = 27r/(i? - 2a), 2j = In A/(i? - 2a) 
is the width of the resonances with the eigenvalue of the 
instability matrix A = [R — a + yjR{R— 2a)]/a. 27 is 
also the classical escape rate 70 of the system. Thus the 
autocorrelation oscillates exactly with period As. Excel- 
lent agreement is found between experiment and theory 
for this 2-disk case. 

In the full space of the 4-disk system, the average 
length of the periodic orbits per period can be esti- 
mated as the average length of the eight periodic or- 
bits, 12, 23, 34, 41, 1234, 1432, 13, 24, where 1,2,3,4 
are the labels of the four disks [^. The mean separa- 
tion between the resonances is approximately given by 
Ai = 27r/[(2 -h 1/V2)i? - (3 + ^/2)a]. For the 4-disk sys- 
tem with one-eighth of the phase space, the mean sepa- 
ration is As = 67r/[(3 + V^)R - 2(2 + \/2)a]. The au- 
tocorrelation will oscillate with approximate period As, 
which indicates the deviation from the semiclassical the- 
ory because of the presence of the periodic orbits in the 
system. Thus the short time behavior is system specific. 
The value of As = 0.35 cm~^ is in good agreement with 
the scale of oscillations in Fig. 2. Future work will focus 
on the interplay between the universal chaotic features 
and the non-universal PO contributions. 

It is interesting to examine the variation of the results 
with increasing n. For n = 00 one obtains the Lorcntz 
scatterer. For n — 20, we find that 7/^4 '-^ 0.05, indicating 
that as n increases the system approaches a closed system 
and the escape rate becomes very small. 

To our knowledge this work represents the first de- 
termination of the wave-vector autocorrelation for an 
experimental system. Our results thus nicely comple- 
ment measurements in semiconductor microstructures 
||12|,p|Jlij| , where however the wave- vector correlation is 
difficult to extract but instead the magnetic field corre- 
lations of conductivity g{A.B) are analyzed. As noted 
previously, the present experiment exactly corresponds 
to an ideal non- interacting electron in a quantum dot, 
with tunneling contacts. Thus Fig.l can be viewed as 
equivalent to conductivity fiuctuations in a quantum dot. 

The present work has provided the first experimental 
realization of the n-disk open billiards problem. The ex- 
periments clearly demonstrate the quantum- classical cor- 
respondence, in that the quantum properties can he ob- 
tained from classical quantities and vice-versa. One of 
the powerful features of the present experiments is the 
(almost unlimited) ability to vary geometry, as we have 
already demonstrated by studying numerous configura- 
tions. Hence there is enormous potential for addressing a 
variety of issues in the quantum-classical correspondence 
problem. In addition open billiards are also a paradigm 
for dissipative quantum systems, a problem of wide im- 
portance. Besides their relevance to atomic and chemical 
physics, the present results also have applicability to sit- 



uations where wave mechanics plays a role, such as in 
electromagnetism and acoustics. 
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